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Abstract 

We present some recent results concerning the persistence of dynamical local- 
ization for disordered systems of n particles under weak interactions. 

1 Introduction 

After more than half a century, Anderson localization continues to attract the interest 
of a broad spectrum of researchers ranging from experimentalists, who finds its effects 
in systems of cold atoms and in photonic crystals [7], to mathematical physicists. Con- 
siderable progress was made in the rigorous methods for the study of the localization 
effects of disorder in the context of the one-particle theory fTTl [141 l6l . More recently 
attention has centered on the role of interactions, and questions related to the persis- 
tence of the localization picture in the presence of inter-particle interactions. We shall 
report here on some progress which was made in that area. 

1.1 The one-particle theory 

The discussion of weakly interacting particles often starts from the approximation in 
which the interactions are ignored. That is, one first considers systems of Fermions, 
or Bosons possibly with on-site repulsion, subject to a fixed potential which includes 
random terms. Such system can be understood in terms of the one-particle theory. 
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For particles moving on a lattice Z d , the single-particle Hamiltonian may take the 
form 

H^(uj) := —A + A V(x; u) , in£ 2 (Z d ), (1.1) 

where —A is the lattice Laplacian and the random potential incorporates a disorder 
parameter A > 0. For the convenience of presentation, all mathematical statements 
in this note are made under the assumption that the random potential V(x; u>) takes 
independent and indentically distributed values on the lattice sites x 6 Z d , with a 
distribution ¥(V X € dv) — g(v) dv of abounded, compactly supported density. 

When talking about localization for operators such as H^(lo), different notions 
have been used and established: 

Spectral localization in an energy regime / c K, refers to the statement that within 
the specified energy regime i?W(cj) has only pure point spectrum, with expo- 
nentially localized eigenfunctions, cf. ifTTl and references therein. 

Dynamical localization in I C K, refers to the statement that any initially localized 
state, which is a wave packet with energies in / c K, will remain exponentially 
localized indefinitely under the time evolution. A convenient sufficient condition 
is that for some A, £ £ (0, oo) (both depending on /, A), 



(DL) E 



where S x <G £ 2 (Z d ) is the S state at x and Pj(H^) is a spectral projection. 

While spectral localization provides some coarse information on the quantum time evo- 
lution via the RAGE theorem, dynamical localization in the sense of (DL) is a some- 
what stronger notion. Historically, spectral localization was the first to be rigorously 
established by building, for d > 1 at extreme energies or large disorder, on the multi- 
scale analysis of Frohlich and Spencer fl2l : cf. ifTTlfPfl . The first proof of dynamical 
localization [2 J relied on the fractional moment analysis of (T). 

It is one of the features of the single particle model that in any dimension 
there is an extreme disorder regime, with A exceeding some finite Ai, in which one has 
complete dynamical localization, meaning that all states are localized and dynamical 
localization (DL) holds with I = R. 

For the single particle model (JTTTJ in d = 1 complete dynamical localization with 
critical disorder strength Ai = was established already in |[T5l . For higher dimen- 
sions, d > 1, an explicit bound on Ai can be found in 0. 



1.2 The influence of interaction on localization: a challenge 

It is not hard to see that the above localization properties are inherited by a system of 
non-interacting particles with a one-particle Hamiltonian corresponding to ( |l.l| i. An 
important question however, is whether such behavior will persist under the addition 
of interparticle interactions. Especially interesting is the situation where there are n 
fermions in a region of volume |A|, with |A| — > oo and n/\A\ — > p > 0. 
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A rather strong claim is being advocated by Basko, Aleiner and Altshuler [8 1, who 
argue that under strong disorder weak interactions do not change the qualitative picture 
of localization as it is seen in the non-interacting model. 

The proposal is rather startling. It includes the claim that if the system of particles 
of an overall positive density is started at a initial state at which the distribution in space 
of particles and energy is far from uniform, its irregularity will persist indefinitely under 
the time evolution. This seems to run against the vague equidistribution principle, 
by which one expects that except under unusual circumstances, such as in the non- 
interacting integrable model, the initial state will evolve in time (in a weak enough 
sense) towards states which maximize the entropy, in a coarse-grained sense, subject 
to the given energy and particle number constraints. 

Persistence of localization is not undisputed among physicist. In particular, long- 
range repulsive interactions are conjectured to have a delocalizing effect lfT6ll . 

Rigorous methods are still far from allowing one to decide whether complete local- 
ization will persist or perish in the presence of interactions. Furthermore, the analysis 
of even a fixed number of particles with short range interactions in the infinite-volume 
limit has presented difficulties, and we will report on some recent progress ifTUl [TT1 [5] 
made in this direction. 

As an side we note that the dynamics of multiparticle systems bear some relation 
to non-linear evolutions. Recent results on that topic are discussed in [ 18 1. 

2 Dynamical localization for multi-particle systems 
2.1 The n -particle Hamiltonian 

We will be concerned with a system of finitely many interacting particles in the random 
potential described above. The n-particle Hamiltonian is given by 

n 

:=5^[-A,- + XV{x j \oj)] +U(x;a), in £ 2 (Z d ) n (2.1) 

acting in the Hilbert space over all configurations x = (x\, . . . , x n ) € (Z d )™. Here the 
last term is a p-site interaction of range £jj < oo: 

p 

U(x; a.) := 

k=l AcZ d :\A\=k 
dia,YnA<£u 

It is described in term of an interaction parameter a = (ai, . . . ,a p ) € W and a 
function U a '■ N'' 4 ' — > K which is bounded by some n-dependent constant, \\Ua\\ < 
c„ < oo, and translational invariant. The latter depends on the number N u (x) :— 
T^j=i 3u,xj °f particles of the configuration x = (xi , . . . , x n ), which are at sites u G A 
in the pattern A c Z d . 

Simply stated examples which are already of interest are short range pair interac- 
tions for which p = 2 and U {u ^ v} (iV„(x), N v (x)) = N u (x)N v (x) 5| u _^| iX . 
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2.2 The main result 



Our main result establishes the existence, for any dimension and any number of parti- 
cles, of a regime in the space of basic parameters of the model, (A, a) G R + x MP, 
for which complete dynamical localization with a uniform localization length occurs 
for up to n particles. The proof gives an inductive algorithm for the construction of 
such localization regimes, albeit at what may be a possibly far too restrictive manner. 

The first proof of spectral localization for two particles was presented in [11], using 
the multiscale approach. The result presented below was derived by different means, 
which allow also rather simple control of the dynamical localization. 

The localization regime covered by the result presented below may be best de- 
scribed in terms of the extreme sets it includes, namely: 

Strong disorder for each aeR p there is A(a) such that the localization regime 

includes the cone in the parameter space R + x M. p where the interaction strengths 
are dominated by a, and the disorder strength exceeds \[ct). 

Weak interactions for any A > Ai, i.e. disorder strength at which the one-particle 
Hamiltonian exhibits complete localization, there are ctj(A) > , j = {1, 

such that effl includes all (A, a') for which \a'j\ < |aj(A)| componentwise. 

Thus, the weak localization region includes some neighborhood of the entire localiza- 
tion regime ( Ai , oo) x {0} of the n particle unperturbed system. In case d = 1, Ai = 0, 
so this includes all positive values of A. 

Having introduced these notions we may present the result of 151 : 

Theorem 2.1 J5J). For each n,p G N there is an open set £„ C R+ x R p which 
includes regimes of strong disorder and weak interactions, for which at some A, £ < oo 
and all (A, a) G £n \ k G {1, . . . , n}, and all x, y G (Z d ) k : 



E 



sup 

.ten 



< Ae~ distH(x ' y)/< ; , (2.2) 



where the exponential decay is in terms of the Hausdorff pseudo-distance between the 
configurations x — (xx, . . . ,Xk) andy = (jjx, . . . , y^) 



dist«(x,y) := max < max dist(:Ej, y), max dist(j/i, x) > (2.3) 

I l<i<k l<i<k J 

As was mentioned, spectral localization follows from (|2.2|i by the Wiener criterion. 



2.3 Comments on the result 



Theorem 2.1 is formulated for distinguishable particles. For Fermions, or Bosons, the 
relevant Hilbert spaces are subspaces of the space which this theorem covers, and thus 
the results apply by restriction. Furthermore, hard core interactions can also be added 
without requiring any modification in the proof of the extended statement. 
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It is worth pointing out that for systems of n > 2 particles some subtleties show 
up in the decay rate seen in ( |2.2[ ). A natural decay rate for systems with permutation 
symmetry is the symmetrized distance: dists(x, y) := min 7r6 s n Y^j=i \ x j ~ V^jl > 
where S n is the permutation group of {1, ...,n}. In contrast, dist^(x,y) is not a 
metric, if n > 2. The Hausdorff pseudo-distance bounds ( |2.2| i allow for the possibility 
that if the initial configuration had some particles within the localization distance from 
each other then such 'excess charge' may transfer among the occupied regions. 




Figure 1: A transition between two configurations which are close in the Hausdorff pseudo- 
distance, but not otherwise. 



One could wonder whether Theorem 2.1 does not immediately follow from the 
existing localization results. After all, a configuration of n particles in % d may be 
regarded as a single particle in % dn evolving under a Hamiltonian which is of the 
form i?(")(u>) = -A(" d ) + U(x) + A V(x; u) with a random potential V(x; io) = 
Y^j=i V{xj\uS). However, the random potential V(x;ui) does not take independent 
values at configurations with one coinciding particle and in this sense is of infinite 
range. The number of its degrees of freedom is only a fractional power (^) of the 
number of configuration. The technical problem which this causes is behind the above 
mentioned limitation of the bound d2.2b. 



3 Outline of the proof 

The proof of Theorem |2. 1 | cannot be fitted in this short summary, but we may comment 
on the flow of the argument [5 1. It consists of three steps: 

1 . Proof of the finiteness of fractional moments of the Green function. 

2. Elucidation of the relation of a multiparticle eigenfunction correlator with the 
Green function's fractional moments. 

3. Inductive construction of domains of uniform localization in parameter space. 

The first is a sine qua non condition for the fractional moment technique for localiza- 
tion JTJ. It is related to the celebrated Wegner estimate, which however does not play a 
direct role in the analysis. The proof is rather easy, as is also the case for the Wegner 
estimate for many-particle systems. 

The second step follows the strategy by which dynamical localization is established 
for single particles Q, for which some adjustments to the multi-particle setup are re- 
quired. 
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Once the first two preparatory steps are taken, the third step constitutes the core of 
the argument. 



3.1 Finiteness of fractional moments of the Green function 

An essential tool for the analysis is the finite-volume, A C Z d , Green function 

G^( X ,y;z):=(6 x ,(H^-z)- 1 6 y ), (3.1) 

with denoting the restriction of ( |2.1| i to I 2 ( A) n . The first step towards exponential 
bounds is to prove finiteness of the fractional moments. In fact, even a conditional 
average makes them finite. 

Lemma 3.1. For any s S (0, 1) there is C < oo such that for any A C 1 d , any two 
(not necessarily distinct) sites U\, £ A and any pair of configurations x,y, which 
have a particle at u\, U2 respectively, 



IS 



(|Gi n) (x,y;z)| S |{^)}^ {ui ,„ 2} ) 



< 



C 

W 



(3.2) 



for all z G C, and A ^ 0. 



Sketch of proof. In its dependence on V{u\) and V(u2), the Hamiltonian is of the 
form: = A+X V{u x ) N Ul +X V(u 2 ) N U2 , where (N u ^) (x) := 6 *w ^( x ) 

stands for the number operator. The assertion is therefore implied by the weak-L 1 es- 
timate from [4]: 



lVN<f>,(ZN + r)M-K) 1 \/MV) 



with operators N, M > and K dissipative. 



> t 



g(0 g( v )d^d V <^-U\\\m, 



□ 



The finiteness ( |3.2) > is the analogue of a Wegner estimate in the multi-scale method. 
As has been noticed in J5J [T3l . it follows similar to the one-particle case from the 
monotonicity used in the above proof. 



3.2 Eigenfunction correlator and its relation to Green function 

(n) 

For the finite-volume Hamiltonian H A , we define the eigenfunction correlator, asso- 
ciated with an energy regime J C t, as the sum 



(x,y;7):= £ | (S x , P {Ey (H™) S y 



(3.3) 



)ni 



where P{e} is a spectral projection. Ignoring some subtleties related to the passage to 
the infinite-volume limit A | Z d (which are discussed in detail in [2, 5 1), the elementary 



bound 



-itH 



(«)\ 



< Q^^x, y; R) shows that the eigenfunction correlator is 
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an important tool for establishing dynamical localization ifTBI Bl. 



The following key lemma states the equivalence of exponential decay in the Haus- 
dorff pseudo-distance (though not in the symmetric distance) of the eigenfunction cor- 
relator and of the Green function's fractional moments. 



Lemma 3.2. The following statements are equivalent: 
1. There is A, £ < oo such that for all x, y £ (Z d ) fe ; 



sup sup E 

/ClAcZ d 



Qi fc) (x,y;I) 



< Ae~ distn{x - y)/i . (3.4) 



2. There is A, £ < oo and s £ (0,1) such that for all x,y£ (Z d ) k : 

sup sup ^ / E [ G^ ] (x,y;E) "1 dE < A e~ dist «( x ^)/« . (3.5) 

I CM AcZ d J I L J 



|/|>1 



=:Ej[...] 



In the proof of the main result we repeatedly change horses between estimates on 
the kernels of G and Q. The kernel Q(x, y) is what is ultimately need, and it is also 
a very convenient tool for the induction step. Yet, G(x,y) has a more convenient 
perturbation theory. 

In the induction step described below one puts together two mutually non-interacting 
subsystems and then turns on the interaction between them. Thus, we start with the sys- 
tem partitioned into J,K C {1, . . . , n}, J n K = 0, and a Hamiltonian of the form 

H[ J K) := H { A J) © H { A K) acting on £ 2 (A)I J I © ^(A)!^. Then: 

• The Green function of the composite system is a convolution, 

dE 

7C-. (3.6) 

2m 



Gi J ' K) (x,y;z)= f G%\xj, yj; z - E) Gf \x K , y K ; E) 
Jc 



involving a contour integral encircling the spectrum of fll ■ Moment estimates 
for the combined system are complicated by the unboundedness of G and possi- 
ble correlations between the eigenvalues of the subsystems. 

• On the other hand, the eigenfunction correlator is bounded and we may use: 

Qi J ' K) (x,y;J) < Q A J) (xj,y 7 ;M) Qf \x K , y K ; K) , (3.7) 
by which the exponential decay of the subsystems passes to the joint system. 
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3.3 Inductive construction of domains of localization in parameter 
space 

In view of the equivalence of the exponential decay in the Hausdorff distance of the 
eigenfunction correlator and of the fractional moments of the Green function, it is 
natural to define localization regimes as follows. 

Definition 3.1. An open subset of the parameter space, C C R+ x R p is said to be a 
domain of uniform n-particle localization if for some s G (0, 1) there exists £ < oo 
and A < oo such that 



(UL) 



1 

SUP SUp t-7 

ICR AcZ d I 1 I 



E 



Gi fc W;£) 



dE < Ae 



d\k 



- dist«(x,y)/£ 



holds for all (A, a) G £, all k £ {I, n}, and all x,y G (Z d ) 

Our aim is to inductively construct starting from the domain C± = (Ai , oo) C 
M+ of uniform one-particle localization, whose existence is guaranteed in Q (and |[T5l 
for d = 1). For the induction step we assume that is a domain of uniform (n — 1)- 
particle localization and proceed by distinguishing two cases: 

1 . Localization for non-clustered configurations. Here we deal with two config- 
uration x, y of which at least one is of diameter (diam(x) := max_y f. \xj—Xf.\) 
comparable with their Haussdorff distance. A key lemma provides the bound: 



sup sup Ej 

ICR ACZ d 

1 



Gi n) (x,y) 



< A cxp ( — — min <j dist^(x, y) 



for all (A, a) G C^-i' (which, strictly speaking, is only valid under the addi- 
tional assumption of being sub-conical, cf. 0). 



max{diam(x) , diam(y) } 



(3.8) 



Figure 2: The solid line is the full and the dot- 
ted lines represent the Green function of the par- 
tially non-interacting system which, by assump- 
tion and due to arguments based on decay 
exponentially. Due to finite range, the interac- 
tion is localized to the shaded strip of width 2£u 
about the diagonal. 



The proof of ( |3.8| > proceeds by breaking the configuration in two parts and removing 
the /nfercluster interaction. Standard perturbation theory then yields an expansion in 
U which in case n = 2, d = 1 is pictorially explained in configuration space in Fig [2] 
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The bound ( |3.8| > is rather crude in its dependence on the particle number. It is the 
main reason why our bounds on the localzation length degrade rapidly with n. 



(ii) Localization for clustered configurations. This refers to configurations with 
diameter less than half their separation. The issue which is to be addressed is 
the possible formation of a quasi-particle which is not constrained by previous 
localization bounds. 

A convenient quantity to monitor is 



flW(L) := sup sup |<9A L | 

m>i 



E 



E 



y£dA L x6C (")(q ; o) 



(3.9) 



where (A; x) denotes the collection of configurations of diameter less than 
L/2 and at least one particle at x. 




• one has at least one particle at 0, and 

• the other one has at least one particle at the 
boundary <9Al of the box; cf. illustration in 
real space. 



Localization bounds for clustered configurations proceed through reseating inequal- 
ities for Bi"\L). The key statement is: 

Lemma 3.3. There exists s € (0, 1), a, A,p < oo, and v > such that 

B^\2L) < j^BW(L) 2 + AL 2l3 e- 2uL (3.10) 

for all (A, a) £ £-^n-v (which is again assumed to be sub-conical, cf. /|5^j. 

Part of the proof of this lemma resembles the proof of finite-volume criterial for 
one particle localization in [3]. In view of the infinite-range correlations of the random 
potential in configuration space discussed at the beginning of this section, a crucial ob- 
servation is the following. When placing two boxes of length L/2 about configurations 
x and y in the above picture, the random variables associated with those boxes are 
independent due to the spatial separation. Moreover, the error when restricting the sum 
in the definition of Bi n \2L) to configurations x,y which have a diameter less than 
L/4 can be controlled by ( |3.8| > yielding the second term on the right side in ( |3. 1Q[ >. 

It is not hard to see that rescaling inequalities such as ( |3. 10| > imply exponential 
decay provided the quantity (L) is small on some scale L. This is the require- 

(v) 

ment which determines £„ . Namely, 
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• in case of strong disorder and arbitrary value of the interaction a. we choose L 
fixed and A small enough. 

• in case of weak interaction and arbitrary value of the disorder A > Ai we appeal 
to ifflE) which im ply that B^\l) -> as L -> oo for a = 0. Since the 
finite-volume quantity Bg (L) is continuous in some neighborhood of ex = 
we choose L and a accordingly. 



4 Some remaining challenges 



While Theorem 2.1 is formulated for interactions of finite range, its proof allows for 
extension to interactions of with exponential falloff. However, it it does not address 
questions about the effects of Coulomb interactions. 

An important outstanding challenge it to resolve the question which is commented 



upon in Section 2.3 does localization persists (with uniform bounds) even for large 
systems at a positive density of particles? 
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